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By definition, tlie Kraus representation of a liarmonic oscillator suffering from tlie environment 
effect, modeled as the amplitude damping or the phase damping, is directly given by a simple 
operator algebra solution. As examples and applications, we first give a Kraus representation of 
a single qubit whose computational basis states are defined as bosonic vacuum and single particle 
number states. We further discuss the environment effect on qubits whose computational basis states 
are defined as the bosonic odd and even coherent states. The environment effects on entangled qubits 
defined by two different kinds of computational basis are compared with the use of fidelity. 
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I. INTRODUCTION 

Since the fast quantum algorithms were developed, 
the quantum information and computation have aroused 
great enthusiasm among physicists because of their po- 
tential applications. The basic unit of the classical infor- 
mation is a bit, which can only have a choice between two 
possible classical states denoted by or 1. The informa- 
tion carriers in the quantum communication and com- 
putation are called quantum bits or qubits. The main 
difference between the qubit and the classical bit is that 
the qubit is defined as a quantum superposition of two 
orthogonal quantum states usually written as |0) and |1) 
called as the computational basis states. It means that 
the qubit has an infinite number of choices of the states 
a|0) + (3\l) with the condition |ap -f = 1. It is ob- 
vious that the superposition principle allows application 
of quantum states to simultaneously represent many dif- 
ferent numbers, this is called the quantum parallelism, 
which enables the quantum computation to solve some 
problems, such as factorization, intractable on a classical 
computer. 

However, a pure quantum superposition state is very 
fragile. A quantum system cannot be isolated from the 
environment, it is always open and interacts with the 
uncontrollable environment. This unwanted interaction 
induces entanglement between the quantum system and 
the environment such that the pure superposition state 
is destroyed, which results in an inevitable noise in the 
quantum computation and information processing. Un- 
like a closed system, whose final state p' can be obtained 
by a unitary transformation U of the initial state p as 
p' ~ U pU\ the final state of an open system cannot 
be described by a unitary transformation of the initial 
state. Quantum operation formalism is usually used to 
describe the behavior of an open quantum system in 
which the final state p' and the initial state p can be 



related by a quantum operation E &s. p' — E{p). In gen- 
eral, the quantum operation £ on the state p can be de- 
scribed by the Kraus operator-sum formalism P, 01 as 
^{p) = ^iJ-P^ji where the operation elements Efj,{Ejj) 
satisfy the completeness relation Ej^E^ = I. Opera- 
tors Ef^ act on the Hilbert space Hs of the system. They 
can be expressed as i?^ = (e^|J7t|eo) by the total unitary 
operator Ut of the system and the environment with an 
orthonormal basis |e^) for the Hilbert space He and the 
initial state ps = |eo)(eo| of the environment. This ele- 
gant representation is extensively applied to describe the 
quantum information processing [3|, |^ . 

In practice, the environment effect on the quantum sys- 
tem is a more complicated problem. There are two ideal 
models of noise called the amplitude and phase damp- 
ing, which can capture many important features of the 
noise "5] . They are applied in many concrete discussions 
to model noise of the quantum information processing, 
e.g., reference j^. As a system, a single-mode harmonic 
oscillator is the simplest and most ideal model to rep- 
resent a single mode light field, vibration phonon mode, 
or excitonic wave. For convenience, we refer this har- 
monic oscillator to a single-mode light field in this pa- 
per. A study of a single harmonic oscillator suffering 
from damping is expected to give us an easier grasp of 
the nature of damping. The Kraus representation of a 
harmonic oscillator suffering from the above two kinds of 
noise has been given by Chuang et al., simply modeling 
environment as a single mode oscillator [3, Q ■ However 
generally speaking, the environment is usually described 
as a system of multimode oscillators for both the phase 
damping and amplitude damping. Milburn et al. have 
also given the Kraus representation of a harmonic oscil- 
lator suffering from the amplitude damping by modeling 
the environment as a system of the multimode oscillators 
using the quantum measurement theory @. Although 
there are many studies about Kraus representation, e.g. 
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in Refs. to the best of our knowledge, there is no 

proof of the Kraus representation for a harmonic oscilla- 
tor suffering from amplitude damping or phase damping 
by directly using its definition when the environment is 
modeled as a system of multiniode oscillators. In view of 
the importance of the Kraus representation in the quan- 
tum information, we would like to revisit this question 
in this paper. We will extend the proof of Chuang et 
al., that is, we will give the Kraus representation using 
a simple operator algebra method by modeling the sys- 
tem as a single mode oscillator and the environment as 
multimode oscillators according to the definition of the 
Kraus representation. 

Our paper is organized as follows. In Sec. II, we give 
the Kraus representation for the amplitude damping case, 
then we discuss the effect of the environment on a qubit 
whose computational basis is defined by odd and even 
coherent states with decay of at most one particle, we 
compare this result with that of a qubit defined by vac- 
uum and single photon number states. In Sec. Ill, we 
give the Kraus representation for the phase damping case, 
and explain its effect. Similar to the analysis in Sec. II, 
we also compare the effect of the phase damping on the 
entangled qubit defined by different states. Finally, we 
give our conclusions in Sec. IV. 

II. AMPLITUDE DAMPING 

A. Kraus representation 

One of the most important reasons for the quantum 
state change is the energy dissipation of the system in- 
duced by the environment. This energy dissipation can 
be characterized by an amplitude damping model. Ide- 
ally we can model the system, for example, a single mode 
cavity field, as a harmonic oscillator, which is coupled to 
the environment modeled as multimode oscillators. Then 
the whole Hamiltonian of the system and the environ- 
ment can be written as 

H = hwh% + h'^uJib\b, + h'^g,{b'^b, + (1) 

i i 

under the rotating wave approximation, where b{b^) is 
the annihilation (creation) operator of the harmonic os- 
cillator of frequency oj, and is the annihilation (cre- 
ation) operator of the ith mode of the environment with 
frequency tUi, Qi is coupling constant between the system 
and the ith mode of the environment. We assume that 
the initial state of the whole system is 

p(0)=P5(0)®P£(0), (2) 

where ps(0) and pe{^) denote the initial states of the 
system and environment respectively. With the time evo- 
lution, the initial state evolves into 

p{t)^U{t)ps{0)®PE{m\t) (3) 



with U = exp{—iHt/h). Since we are only interested in 
the time evolution of the system, we perform a partial 
trace over the environment and find the reduced density 
operator of the system as 

Ps{t) - Tr£ {U{t)ps{0) ® PE{0)UHt)} (4) 
^Y.(^h}\U{t)psiO)®PEmHt)\{h}) 

where = \ki ■ ■ ■ ki ■ ■ ■ ) = 11^(8) |fci)} is an orthonor- 

mal basis of the Hilbert space He for the environment, 
and ki denotes that there are ki bosonic particles in the 
ith mode. We assume that the environment is at zero 
temperature, so its initial state is in the vacuum state 
|{0}) = Ili (g) \0i) = I • • • • • • ) of multimode harmonic 
oscillators. The Hilbert space He can be decomposed 
into a direct sum of many orthogonal subspaces Hk as 
He = ®feLo with an orthonormal basis 
which satisfies the condition ki = k for each sub- 

space Hk. So we can regroup the summation in Eq.l^J 
and define the quantum operation element as follows 

k 

Au{t) = Y,' {{h}\u(t)m) (5) 

which acts on the Hilbert space Hs of the system. Here- 
after, ^ ' stands for summation under the condition 
fci — k; Ak{t) means that there are k bosonic par- 
ticles absorbed by the environment through the evolu- 
tion over a finite time t, and Aj,{t) denotes the Hermi- 
tian conjugate of Ak{t). Then Eq.Q) can be rewritten 
by virtue of the operation elements Ak {t) and its Hermi- 
tian conjugate Al^t) as a concise form called the Kraus 
representation 

oo 

ps{t)^Y.Mt)psmi{t). (6) 

k=Q 

It is very easy to check that the elements Ak{t) satisfy 
the relation Al.{t)Ak{t) = I. In order to further give 
the solution of the operator Ak{t), we first use the num- 
ber state |n) of the system to define an orthonormal basis 
{\n)} of the Hilbert space Hs of the system. Then the 
matrix representation of the operator Ak {t) can be writ- 
ten as 

Ak{t)^Y.^rn,nit)H{n\ (7) 

with the matrix element 
fe 

Atjt) = Y,' {m\{{h}\U{t)m)\n). (8) 

{k.} 

Because the state [/(t)|{0})|rt) can be written as 
C/(t)(6t)"/V^|{0})|0) = pi-t)]"/V^.m)\0). Using 
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the Wigner-Weisskopf approximation [9|, the operator 
(t) can be obtained by solving the Heisenberg operator 
equation of motion for Hamiltonian Q as 

oo 

b\t)=u{tp{Q) + Y,v,{t)bl{0), (9) 

1=1 

where u{t) = exp (— + iwt) and the Lamb shift has 
been neglected. The damping rate 7 is defined as 7 = 
27re(ci;)|(7(w)p with a spectrum density e{uj) of the envi- 
ronment. The details of the time dependent parameter 
Vi (t) can be found in 0, ^ . We can use Eq. to ex- 
pand U{t)\{0})\n) as many terms characterized by the 
particle number k lost from the system, and each term 
with fixed k can be expressed as 

[.nor- ^ t i'jj^. 

k I 

^\{{<SMS^))^ m)\n-k). (10) 

9=1 

By using Eas.lf7|)- (fTn|l . we can obtain the product of the 
matrix elements of „ and its Hermitian conjugate as 

^™,j^'™,„)^ ^[iym^''-''^[^-mTK..n-k (11) 

with ri{t) — e""*"*, where we have used the completeness 
of the subspace with fixed k 

k 

E'i{fc*})«^^}| = i' (12) 

{fc,} 

and condition |w(OP + Sj — which comes from 

the commutation relation [&(t), 6^^(t)] — 1. It is clear 
that the non-zero matrix element A'^ „ (t) must satisfy 
the relation n — m = fc, and all contributions to the 
elements of the quantum operation Ak{t) come from the 
states \n) satisfying the condition n > k. Then we can 
take each matrix element A'^ „(t) as a real number and 
finally obtain HQ 

Mt) = f;y^[ry(t)]("-^-)/2[l-rKt)]'=/2|n-fc)(n|. 

(13) 

It is not difficult to prove that J2k M:i^)^kit) — 1. In 
Eg. lfT^ . one finds that(l — ?7(t))'''/^ is the probability that 
the system loses k particles up to time t, or the prob- 
ability that the state \n) is undecayed corresponds to 
[77(i)]^"~'^^/^ for k particle decay process up to time t. 
For the convenience, we will write r](t) as 77 in the follow- 
ing expressions. 

Using a single qubit p defined by the bosonic number 
state {|0), |1)} as the computational basis states, we can 
simply demonstrate how to give its Kraus representation 



when it suffers from the amplitude damping. By Ea. (|13|l . 
it is very easy to find that the quantum operation on 
qubit p includes only two operational elements Aq and 
Ai, that is 

£{p)^AopAl + A,pA\ (14) 

with 

Ao{t) - |0)(0| + V^|1)(1|, (15a) 
A,{t) = yT^|0)(l|, (15b) 

where y/1 — rj is the probability that the system lose one 
particle up to time t. 

B. Effect of amplitude damping on qubits 

In this subsection, we will demonstrate an application 
of the above conclusions about the Kraus representation. 
We know that two kinds of logical qubits, whose compu- 
tational basis states are defined by using the bosonic even 
and odd coherent states or vacuum and bosonic sin- 
gle particle number states, are accessible in experiments. 
For brevity, we refer to qubits defined via the even and 
odd coherent states as the (Schrodinger) cat-state qubits 
by contrast to the Fock-state qubits defined by vacuum 
and single-photon states in the following expressions. It 
has been shown that the bit flip errors caused by a single 
decay event, which results from the spontaneous emis- 
sions, are more easily corrected by a standard error cor- 
rection circuit for the logical qubit defined by the bosonic 
even and odd coherent states [l^ than that defined by 
the bosonic vacuum and single particle number states. 
We know that the entangled qubit takes an important 
role in the quantum information processing. So in the 
following, we will study the effect of the environment on 
entangled qubits defined by the above two kinds of com- 
putational basis states when they are subject to at most 
one decay event caused by the amplitude damping. Let 
I a) be a bosonic coherent state (we denote this single 
bosonic mode as a single-mode light field, but it can be 
generalized to any bosonic mode, for example, excitonic 
mode, vibrational mode of trap ions and so on). In gen- 
eral, if a single-mode light field, which is initially in the 
coherent state |a), suffers from amplitude damping, the 
Kraus operation element H13|) changes it as follows 

Ak\a)^e 2 -= \y/va), (16) 

and the normalized state, which is denoted by A/'[Afe|a)], 
can be written as 

^[Ak\a)] = IV^a), (17) 

which means that the coherent state remains coherent 
under the amplitude damping, but its amplitude a is 
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reduced to ^Jr\a due to the interaction with the environ- 
ment. 

Now we define a logical zero-qubit state |0) l and a log- 
ical one-qubit state |1) l using the even and odd coherent 
states as 



|0)l = iV+(|a) + |-a)) 
= 7V_(|a)-|-a)) 



(18) 



with N± = (2 ± 2e~^l"l^)~^ and subscript L denotes 
the logical state. If we know that the system loses at 
most one photon with the time evolution but we do not 
have the detailed information on the lost photons. Then, 
we only need to calculate two Kraus operators Aq and 
A\ corresponding to no-photon and single-photon decay 
events, respectively. 

No-photon decay event changes the logical qubits |0)l 
and as follows 
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AA[Ao|0)l] =7V;(IV^a) + |-0^a)), (19) 
AA[Ao|l>L]-7V:(IV^a)-|-V^a», (20) 



where |0) and |1) are normalized states of Ao|0)i 
and Ao|1)l with the normalized constants N'j^ — 

{2 ± 2 exp(— 277|ap)} . It is obvious that no-photon 
decay event reduces only the intensity of the logical qubit, 
leaving the even and odd properties of the logical qubit 
states unchanged. A simple analysis of how the noisy 
channel affects the original quantum state can be made 
by calculating the fidelity defined as the average value 
of the final reduced density matrix p(f) with initial pure 
state IV'), that is /(i) = (V'IpI'/')- The fidelity /^(t) of 
the qubit state H18() with no-photon decay event can be 
obtained as 



(eV^|a|^ _(_ (_l)ie-v^|a|^)2 



(el"l' (-l)*e-l"l')(e"l"l' + (-l)''e-''l"l') 



with i = 0, 1, which denotes zero-qubit state or one- 
qubit state. But for a single-photon decay event, we can 
have the following 



|l)=AA[Ai|0)i], 



|0)=AA[^i|l)i] 



(22) 



We find that the single photon decay event fiips an even 
coherent state to an odd coherent state with the reduc- 
tion of the amplitude and vice versa. It is very easy to 
prove that 

|T) = AA[A(2„,+i)|0)l] =AA[A2„|1)l], (23) 
|0) = AA[^(2™+i)|1)l] =AA[A2„|0)l] (24) 

with m = 0, 1, 2 • • • , which means that the even number 
of photons decay event only reduces the intensity of the 
logical signal, and keeps the even and odd properties of 
the logical state unchanged; but the odd number of pho- 
tons decay event flips a qubit causing error and reduces 
the intensity of the qubit signal. If the coherent ampli- 
tude of the cat states is infinitely large, that is, |a| oo. 



then we find that the fidelity fiit) k, 1 with a larger ry, 
so under this condition, we can say that the no-photon 
decay event essentially leaves the logical qubit states un- 
changed and the single-photon decay event causes a qubit 
flip error. It means that the use of the even and odd 
coherent states as logical qubit states |0)l and |1)l is 
better than the use of the vacuum state |0) and single- 
photon state |1) as logical qubit states in the amplitude 
damping channel with a few photons loss. Because the 
single-photon decay event changes vacuum and single- 
photon states to ^i|0) = and A\\\) = |0), the single- 
photon decay event is an irreversible process for vacuum 
and single-photon states. 

An entangled pair of qubits, whose computational ba- 
sis states are deflned by the even and odd coherent states 
(|18|l . can be written as 

|*) = -i={|0)L|l)L + |l)L|0)4. (25) 

On the basis of the above discussions, we can study the ef- 
fect of environments on entangled qubits when each mode 
loses at most one photon. For simplicity, two modes are 
assumed to suffer from effect of two same independent en- 
vironments. There is no direct interaction between two 
systems. After the environment performs the following 
four measurements {Eq = Ao(>^Ao, Ei — Aq^:^ Ai, E2 ^ 
Ai ^ Aq, E3 = Ai ^ Ai} with Aq an Ai determined by 
Eg. lfT^ . the entangled qubit becomes the following 
mixed state 



Pit) 



J:loE^p{0)E} 



|ii)(ii| + |66)(ii 



= iV<i2p|a|2(— 100)(00| 



N': 



|11)(00| 



-(l+p2|«|4)(|0i) + |i0))((0i| + (i0|)} (26) 



(21) 

with p = 1 — 7] and the normalized constant 



-p\a\ 



1 



-4jj|aP 



(27) 



The fidelity Fi {t) of the entangled qubit I^E") with at most 
one photon decay for each mode can be calculated as 

F,{t) = {^Wm (28) 
(14- |a|V) csch(2|a|2)sh2(2^|a|2) 
" 2|a|2p ch(2|a|2r;) + (1 + |a| V) sh(2|a|2ry) ' 

If each mode can dissipate any number of photons, but 
we know nothing about the details of dissipation, then 
we must sum up all possible environment measurements 
on the system, obtaining the reduced matrix operator of 
the system as 



P'it) 



1 



2 - 2e-4|"l 
M 



.{|/?)|/3)(/3|(/3K|-/?)|-/?)(-/3|(-/3|} 



2 _ 2e-4|"l' 



{\-P)\-p){m + \pm{-P\{-m (29) 



5 



with M = e-'il^l'Ci-'^ "'*) and = \u{t)a). Then we 
obtain the fidehty F2 (t) as 

F2{t) ^ch{2\a\^p)csch^{2\a\^)sh^{2^\a\^). (30) 

We can find that the logical state |0) l and in Ea. H25|) 
can be reduced to the vacuum state |0) and single-photon 
state |1), respectively in the limit of the weak light field 
|a| ^ 0. It is easy to check that 

Ao\0) = |0), ^|1)=V^|1), (31) 
Alio) = 0, A,\l) = yr^lO) (32) 

which means that the |0) qubit state is invariant when 
no-photon decay event happens, but the amplitude of 
the |1) qubit state is reduced to y/rj. If the system leaks 
one photon, then the |0) qubit state vanishes and the |1) 
qubit state returns to the |0) with the probability ^1 — rj. 
The entangled qubit, whose computational basis states 
are defined by vacuum and single-photon states, can be 
written as 

|*) = i=(|0>|l) + |l>|0)) (33) 

which can be obtained by setting \a\ in Ea. (|25|l . 
We can use the same step to obtain the fidelity F{{t) 
corresponding to Ea. (|33|l when both modes are subject 
to the amplitude damping as 

m)-vit), (34) 

where we assume that the damping is the same for the 
two modes. When the weak light field limit \a\ ^ is 
taken, we find that Ea. H28|l and Eg. 1)30(1 can be written 

as 

Flit) = r^ + |(l -577 + 37,2 ^773)|a|4 + 0(|«|8), 

F^it) = rj + ^77(1 - 47; + 3772)|a|4 + Oi\af ). (35) 

It is clear that all the higher order terms of |ap can be 
neglected in the limit |a| 0, so Eq.(|2HIl and Eg. lfSUI) 
can return to F[{t) in the weak field approximation. By 
recalculating the difference 

Flit) - F2it) = -177(1 - rj)'\a\^ + Oi\af), (36) 

it is clear that fidelity Fi(i) is less than F2(t) in the 
low-intensity regime. These predictions are confirmed 
by the exact evolution of fidelities depicted in Fig. 
for \a\ < 1. To compare the fidelities in the short-time 
evolution, we expand the difference of Eqs. H28|l and H30|) 
in power series in Tt as follows 

Fiit)~F2it) - -^|a|6coth(2|an(7t)3 + 0{(T-i)4} 

(37) 

which shows that Fi it) is less than F2 it) for arbitrary in- 
tensities in the short-time regime although the difference 
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FIG. 1; Difference of fidelities: Fi ~ rj (solid curves), given 
by Eqs. 12811 and 1341 . and F2 — r] (dashed curves), given by 
Eqs. liniU and versus scaled time for different values 

of amplitude ct. 

might be very small as presented in Fig. ^ for |a| — 2,4. 
In the long-time limit, we can expand the difference in a 
series of rj resulting in 

Fiit)-F2it) = ^J^[l-H|a|4-2|apcoth(2|ap)]77 
sh(2|Q;|-^) 

+0{ff} (38) 

which is negative only for the intensity |ap < 1.1997 • • • . 
For higher intensity, we find that |1) « and |0) ~ 
\Q)L^ and p — \ — rj will become small with the time evo- 
lution, then the fidelity Fiit) will become greater than 
the fidelity F2it) in some regimes. But in the low inten- 
sity limit we find that |1) « |0) and |0) « |0), then the 
fidelity Fiit) will become less than i^2(i)- 

III. PHASE DAMPING 

A. Kraus representation 

A state can be a superposition of different states, which 
is one of the main characteristics of the quantum mechan- 
ics. The relative phase and amplitude of the superposed 
state determines the properties of the whole state. If the 
relative phases of the superposed states randomly change 
with the time evolution, then the coherence of the quan- 
tum state will be destroyed even if the eigenvalue of the 
quantum system will be changed. This kind of quantum 
noise process is called the phase damping. We can have 
the Hamiltonian of a harmonic oscillator suffering from 
the phase damping as 

H = hwb'^b + cj,6j&, + h^Xtb^Kbt + bl), (39) 

i i 
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where b{b^) is the annihilation (creation) operator of the 
harmonic osciUator with frequency a;; and is the 

annihilation (creation) operator of the ith mode of the 
environment with frequency uji. Xi is a coupling constant 
between the system and the ith mode of the environment. 
We can solve the Heisenberg operator equation of motion, 
and very easily obtain the solution of the system operator 
b{t) as 



b{t) = 6(0) exp —iujt 



bj{t) = bj(0) exp ^—iujjt — ixjtb^b^ ■ 



(41) 



We also assume that the system-environment is initially 
in a product state p{0) = ps{0) Pe{0)- We apply the 
time-dependent unitary operator U{t) = e^'^*/'' with 
the Hamiltonian H determined by Eg. 1)39(1 to the state 
|{0})|to) as follows 



UmO})\m) 



mm = 



^{-tmcjt-i T. X,mt{b] (-*)+&, (-*))} ^^"(0) 



|{0})|0)(42) 



By assuming that the environment scatters off the quan- 
tum system randomly into the states {\ki ■ ■ ■ ki ■ ■ ■)} with 
the total particle number fc, then the Kraus operators can 
be defined as follows 



{k^} 



(43) 



with its Hermitian conjugate -Using Eq (|42|l and Eq 
(|43|l . the reduced density operator of the system can be 
expressed as 



(44) 



fc=0 



with 



F,(t)^Pfc(r) = f;exp|-inV| 



(n^r^) 



fc! 



\n){n\ 



with p = |7/^)('(/'|, and 



Pfe(r) =<5feo|0)(0|+exp 



-^rH^IDdl, (47) 



where 5ko is Kronecker delta. By contrast to the case of 
the amplitude damping in which k can only take and 1, 
in the case of the phase damping, k can take values from 
to oo, which means that the system can be scattered by 
any number of particles in the environment. In this sense, 
we can say that there is an infinite number of quantum 
operation elements Pfe(t) acting on the qubit when 
it suffers from a phase damping. We find that Pk^o(t) 
makes the state |0) equal to zero, so we can rearrange all 
operation elements Pk^o{t) as one group, and redefine 
two Kraus operators as 



En 



exp^-ir2^|l)(l| 



= v/l-cxp{-r2}|l)(l| 



(48) 



which is the form of the reference P). It is evident 
that EqEq + EiE\ — I, and Eq means that the en- 
vironment returns to the ground state after it is scat- 
tered by the system. Unlike the amplitude damping, Ei 
means that the environment is scattered to another state, 
which is orthogonal to its ground state. It does not mean 
that only one particle in the environment is scattered by 
the system. We find that the probability of a photon 
from the system being scattered by the environment is 
^1 - exp{-r2}. 



B. Phase damping effect on qubits 

In this subsection, we will further discuss the phase 
damping effect on a qubit whose computational basis 
states are defined by the bosonic even and odd coher- 
ent states or vacuum and single photon number states. 
If the coherent state \a) of the system is scattered by k 
photons of the environment, then it changes as 



n=0 



(45) 

where t — t^/V is a rescaled interaction time and V = ^ ^ 

Y.j \Xj?- \/l - exp{-n2T2} can be interpreted as the Pfc(r)|a) = ^= ^ exp{--(n^r^ -I- \a\'^)} 
probability that n particles from the system are scattered ' ^' 

by the environment. It is clear that a sum of all Pfc 

satisfies the condition X^felo -^fe(*)-^fc(*) = ^- ^^^^^ 
use qubit as a simple example to investigate the phase 
damping effect on it and give its Kraus representation. 
We assume that our harmonic system is initially in |^) = 
q;|0)-|-/3|1) and suffers from the phase damping. By virtue 
of Eq. (|^ , it is easy to find that the quantum operation 
on the qubit can be expressed as 



(49) 

We are interested in the change of the off-diagonal terms 
for the system state after it scatters the photons of the 
environment, for k photons scattering, we have 



2k 1 

Pk{T)\a){P\Pl{T) = ^— exp{--(|ap -I- W)} (50) 



E 

fc=0 



Pk{T)pPl{T) 



(46) 



E exp 



n,m— 



-{n +m)\ == \n){m\. 
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As the system can scatter an arbitrary number of photons 
in the environment, we need to sum all P]^{t), getting 



Y,Pu{r)\c^){P\Pl{r) = cxp{--(|ap + W)} 



\n) {m\ 



(51) 



In order to get a clear illustration of way how the phase 
damping affects the logical qubit states, we can write out 
reduced density matrix of logical states \0)l and in 
Ea. (fTH|l as follows 



Po(-r) 
Pi{t) 



g-2T (n-m) ^2n^*2m 



\2n){2m\ 



(52) 



-2r^(n-m)^ 2n+l *2m+l 



n.ni 



sh| 



ar y/ {2n + iy.{2m - 



IV. 



|2n+ l)(277i+ 1| 
(53) 



It is obvious that when the qubit states are subject to 
phase damping, although the even and odd properties 
are not changed, all of the off-diagonal elements of the 
reduced matrix pk{t) with (k — 0, 1) tend to zero with 
the time evolution. In contrast to the amplitude damp- 
ing, the qubit states suffering from the phase damping 
are no longer pure states even though no photon is lost. 
Because the relative phases of different superposed ele- 
ments of the coherent state have been destroyed by the 
random scattering process of the environment. The fi- 
delity fk{t) (k = 0, 1) of the logical qubit states (|18|l can 
be obtained as 



h{T) = 



1 



ch^laP - k 



,\i{n+m+k) -2{n-m) 



{2n + k)\{2m + k) 



(54) 



where /^(t) < 1(A: = 0, 1) with the equality sign hold- 
ing only in the limit |a| so that |0)l ~ |0) and 
~ But if we directly choose the vacuum state 
|0) and single photon state |1) as the logical qubit states, 
and they independently go into the phase damping chan- 
nel, the fidelity is kept as one. It is because the phase 
damping only changes the off-diagonal terms. It means 
that in the phase damping channel, the use of the vac- 
uum state |0) and single-photon state |1) as logical qubit 
states is better than the use of the even and odd coherent 
states as logical qubit states |0)l and |1)l. The fidehties 
of the even and odd coherent states suffering from phase 
damping are compared in Fig. We find that in the 
low coherent intensity of |ap, the even coherent state 
can keep a better fidelity than that of the odd coherent 
state. But the fidelities for the even and odd coherent 
states approach to each other with increasing intensity, 
so the states have no difference for the loss of information 
in the high coherent intensity. 

To show this property analytically, we find that one 
summation in Eq. (|54|l can be performed leading us to 




0.5 1 



FIG. 2; Evolution of fidelities /o(t) (solid curves) and /i(r) 
(dashed curves) are given by Eq. 15411 for different values of 
intensity |a[^. It is seen that the curves for /o(r) and /i(t) 
coincide if lap > 4. 



relation {k = 0, 1) 
1 



fk — 



{h{2\a\') + {~lfM2\a\') 



2(ch2|a|2-A:) 

OO 

+2^6-^'^'-' [l2d{2\a\^) + {-lfj2d{2x)]\bb) 



d=l 



where J2d{x) is the Bessel function and hdix) is the 
hyperbolic Bessel function. If |ap ^ 1 (on the scale 
of Fig. 2, for lap > 4) then Id{2x) > Jdi2x) and 
(ch^|ap — k)~^ « e^^l"! , which implies that both fi- 
delities fk are approximately equal to 



/o « /i « 2e-2|"l^ /o(2|c 



2E' 



-2d-' 



l2d{'2\a\') 



(56) 

Assuming also long scaled interaction times, = Tt^ ^ 
1, the fidelities can be approximated by compact formulas 

/o«/i « 2e-2|"l'[/o(2|an + 2e-2r^/i(2|an] 

« -^[l-f2exp(-2r2)]. (57) 



We can also study the effect of the phase damping on 
the entangled qubit H25() based on the above discussions, 
and the relevant calculations are straightforward. Here, 
we only give a fidelity F{t) for two modes subjected to 
the phase damping as 



F{t) 




2(2ri+2m+l) -2(n-m)(n-m+l)7 



(58) 

\ 2 



(2m)!(2n+ 1) 



where /o(t) and /i(t) are given by Eq. H54|) . Assuming 
that both channels have the same phase damping con- 
stant F. Eq. H58|l simplifies to the following single-sum 
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FIG. 3: Evolution of fidelities F'{t) (dashed curve), given by 
Eq. 1621 . and F{t) (solid curves), given by Eq. 11581 . for 
different values of intensity \a\^ ■ 



formula 



(59) 



^g-2d(<i+l)r" 



X I > e--^-^'''l2d+i{2\a\^)+Ii{2\a\^) 



which helps us to find approximate compact-form so- 
lutions. E.g., for either long scaled interaction times, 
T » 1, or low intensities, |ap <^ 1, the fidelity F{t) 
reduces to 



M2lc,l')]\ (60) 



On the other hand, by assuming high intensity, |ap 3> 1, 
together with r ^ 1, one obtains the following approxi- 
mation 



F(r) 



1 



27r|a|' 



[l-f 4exp(-r2)] 



(61) 



If the qubit states are defined by the vacuum state 
|0) as logical zero state and single-photon state |1) as 
logical one state, then for the entangled qubits (|^ . if 
both of them go into phase damping channels with the 
same damping constant F, we can obtain the fidelity as 



(62) 



In Fig.®, we compare the fidelities for the Eas. (|58|) 
(|62|l . We find that their fidelities are almost the same 
when the odd and even coherent qubit states are in the 
low intensity limit. But when the intensity increases, the 
fidelity given in Ea. (|^ is less than that in Ea. (|^ . It 
shows that in case of the phase damping, the odd and 
even coherent states are not good coding states. 



IV. CONCLUSIONS AND DISCUSSIONS 

Using a simple operator algebra solution, we have de- 
rived the Kraus representation of a harmonic oscillator 



suffering from the zero temperature environmental effect, 
which is modeled as the amplitude damping or the phase 
damping. It is worth noting that we could not derive a 
compact form for the Kraus representation when the en- 
vironment is initially in the thermal field, because it be- 
comes very difficult to distinguish the number of photons 
detected by the environment and that absorbed by the 
system due to the existence of the thermal field, which 
can excite the system to higher states fl^ . 

As examples and applications of a system in a zero 
temperature environment, we first give a Kraus represen- 
tation of a single qubit whose computational basis states 
are defined as bosonic vacuum and single particle number 
states. We further discuss the environment effect on the 
qubits whose computational basis states are defined as 
the bosonic odd and even coherent states. We find that 
when the system suffers from the amplitude damping, the 
loss of even number of photons leaves the qubit, whose 
computational basis is defined by the even and odd coher- 
ent states, unchanged, but the loss of an odd number of 
photons changes the even or odd properties of the qubit. 
If the system loses a few photons and the intensity of 
the coherent states is taken infinitely large, then we can 
roughly say that the loss of an even number of photons 
keeps qubit unchanged, but the loss of an odd number 
of photons causes a bit flip error. Such an error can be 
corrected by some unitary operation. But if the compu- 
tational basis is defined by the vacuum and single-photon 
states, a single-photon loss is an irreversible process, we 
cannot find any unitary operation to correct this error 
resulting from a single photon loss. So in the amplitude 
damping channel, the use of the even and odd coherent 
states as logical qubit is more suitable than the use of the 
vacuum and single-photon state as logical states 0, 0| . 

When the qubit states, whose computational basis is 
defined by the even and odd coherent states, enter the 
phase damping channels, the even and odd properties 
will not be changed, however all off-diagonal terms of 
the qubit density matrices gradually vanish. Then the 
originally pure states change into mixtures of states be- 
cause of the random scattering of the environment on the 
system even without loss of photons. It is much more dif- 
ficult to correct such an error. So the use of the even and 
odd coherent states as logical qubit states cannot solve 
the problem of the phase damping. 
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